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Time-dependent close-coupling calculations of the ionization and excitation cross section for an- 
tiproton collisions with molecular hydrogen ions are performed in an impact-energy range from 0.5 
keV to 10 MeV. The Born-Oppenheimer and Franck-Condon approximations as well as the impact 
parameter method are applied in order to describe the target molecule and the collision process. 
It is shown that three perpendicular orientations of the molecular axis with respect to the tra- 
jectory are sufficient to accurately reproduce the ionization cross section calculated by [Sakimoto, 
Phys.Rev. A 71, 062704 (2005)] reducing the numerical effort drastically. The independent-event 
model is employed to approximate the cross section for double ionization and production in 
antiproton collisions with H2. 



PACS numbers; 34.50.Bw 



I. INTRODUCTION 



A remarkable progress in the understanding of in- 
teractions between antiprotons {p) and atoms has been 
achieved over the last decades. Although the number of 
antiproton collision experiments [J, [3, UH, [^, 0, 0] is lim- 
ited due to the required effort for the production of low- 
energy p, a large amount of theoretical studies employing 
a variety of different methods have been performed focus- 
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get s but also other targets like alkali-metal 
|32 | atoms have been considered. Among these are full 
quantum-mechanical treatments for H as (loj and fully- 
correlated two-electron calculations for He, e.g., (20l. |2]|. 
They provide single-ionization cross sections in agree- 
ment with experiment for _E > 40 keV. For lower energies, 
however, some open issues still remain in the comparison 
among the theoretical results and with the experiments. 
In the case of double ionization of He differences have 
been observed between proton {p) and antiproton exper- 
iments. The measured data for p are larger by a factor 
two even for energies around 1 MeV, where the single- 
ionization cross sections for p and p are virtually the 
same, and were reproduced by theoretical calculations 

For p collisions with molecules experimental p -I- H2 
data for ionization cross sections [l|, [33| and for stop- 
ping power [sl, [111 were measured. As for the double 
ionization of helium targets, a considerable difference in 
the production of H+ ions was observed between p and 
p [So , Wf\ impacts. Theoretically, however, only little 
has been investigated for p impacts on molecular targets. 
Very recently, the single-ionization and single-excitation 
cross sections [38)] as well as the stopping power for p 
+ H2 were calculated using a one-electron one-center de- 
scription of the H2 target [40| . The experimental p results 
could only be reproduced for impact energies i? > 90 
keV. The findings suggest that for lower energies molec- 
ular as well as electron-correlation effects become impor- 



tant and have to be considered. 

An elaborate calculation of the ionization process for 
molecular targets was performed recently by Sakimoto 
(4l[ being a pioneering work on p -\- H^ collisions. The 
calculations were performed using a discrete variable rep- 
resentation method (DVR) in which the wavefunction is 
directly calculated on grid points constructed from the 
zeros of orthogonal polynomials. The author examined 
closely the dependence of the ionization cross section on 
the internuclcar distance of both nuclei and the rela- 
tive orientation of the molecular axis with respect to the 
trajectory of the antiproton. Thereby, it was possible 
to present cross sections for ionization independent of a 
fixed orientation of the molecular axis and internuclcar 
separation. However, the author considered the effort for 
these calculations due to the many degrees of freedom as 
being extremely time-consuming especially in the case 
that also excitation cross sections are considered. The 
calculation of the latter was postponed although the em- 
ployed method is in principle capable of describing exci- 
tation. 

It is therefore one aim of the present work to reduce 
the amount of time needed for the computation oip colli- 
sions with molecules. This is not only important for the 
determination of HJ excitation cross sections but even 
more a prerequisite for calculations of more complex sys- 
tems like p -I- H2 which are in the focus of the ongoing 
research. In order to decrease the computational effort 
a number of different actions are taken into account in 
the present work. These are an appropriate basis repre- 
sentation using eigenstates of the unperturbed Hj ion, a 
reduction of the number of different orientations of the 
molecular axis which are calculated, and the use of the 
symmetry of the collision system to reduce the number 
of coupled equations. 

The information on the H+ production in p + H^ col- 
lisions can be used for the interpretation of the H+ pro- 
duction in p -|- H2 collisions which was measured one 
and a half decades ago [l| but is still not understood 
theoretically. In [4l| the idea was discussed to use the 
concept of a two-step sequential ionization process to 
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explain the measured p + Ho data in analogy to what 
was done by Janev et al. [12] and Wehrman et al. 
in order to describe double ionization in p + He colli- 
sions using an effective single-electron descriptions. The 
underlying idea is that in a first step the target is ion- 
ized and one electron is ejected. In a second step the 
projectile interacts with an ionic rather than a neutral 
target reducing the probability for double ionization. It 
turned out that Wehrman et al. were fairly successful 
with this independent-event model (lEV) in reproducing 
the measured antiproton double-ionization cross section 
for He targets by using the product of the single-electron 

probabilities pP^i ^^'^ Pkm ^'^^ ionization of He and He+ , 
respectively. 

In section [H] the present method and the used approxi- 
mations as well as relevant symmetries are described. In 
section Hm the convergence behavior of the results is stud- 
ied. Afterwards the calculated ionization and excitation 
cross sections for p + HJ are presented and discussed. 
The obtained data are used to estimate double ionization 
and H^ production in p collisions with H2 molecules. Fi- 
nally, section Hvl closes with a summary and conclusions. 

Atomic units are used unless stated otherwise. 



II. METHOD 



For the description of the electronic wavefunctions 

i.„ax(n) 

V'jyn(r;i?nuc)- J2 Pm{r;Rnnc)Yi'\u) (3) 
i=i,„i„(M,n) 

a one-center expansion around the midpoint of the inter- 
nuclear axis is chosen, where r and uj are the radial and 
angular variables of the electron, respectively. In the ex- 
pansion in Eq. ^ only even or odd values of / contribute 
depending on whether H is gerade or ungerade, respec- 
tively, and Irain > M . Notc, Imin and Imax are merely 
basis-set parameters. The angular part is described with 
spherical harmonics Yj^'^ (lj). The z'-axis of the molecule- 
fixed space is chosen along the internuclear axis. The 
radial part pfJi{r) is expanded in a B-spline basis of the 
order 8. The radial equation is solved in a finite box with 
a radius of 100 a.u. using fixed boundary conditions by 
what bound as well as discretized continuum states are 
obtained. The appropriate representation of the contin- 
uum is an advantage of the used basis expansion which 
was already successfully employed before for atomic tar- 
gets j31i] . 

The electronic structure code which solves Eq. ([2]) uses 
a one-center approximation of the molecular potential 



^(r,i?nuc) = 
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A. Target description 



which is expanded using the relation found by Legendre 



In the Born-Oppenheimer approximation the total 
wavefunction for the Hj molecule separates into the 
product 

i^ir, R„,,) = ^I2^^Yp{e, <f )V(r; , (1) 

where Xi'j are the eigcnfunctions of the molecular vi- 
bration, YJ^ the spherical harmonics, and {v, j, m) the 
vibrational and rotational quantum numbers. Rnuc = 
(i?nuc, ©5 *!') and r are the position vectors of the nu- 
clei and the electron, respectively. The wave func- 
tion ^{r; Rnuc) satisfies the electronic part of the time- 
independent Schrodinger equation 

He ^j^n(r; ^nuc) = 4n(^nuc) V-ilnlr; -Rnuc) (2) 

for an unperturbed molecule at a fixed internuclear dis- 
tance -Rnuc, where is the sum of the potential and 
the electronic part of the kinetic operator. In contrast 
to atomic targets which are spherical symmetric the 
two-center mono-atomic molecule ion H^obeys different 
molecular symmetries. Instead of the atomic quantum 
numbers the electronic part of the HJ eigenstates tp can 
be characterized by N, M and H, where H is the per- 
mutation symmetry with the values g for gerade and u 
for ungerade symmetry, M is the projection of the an- 
gular momentum on the internuclear axis, and N is the 
principal quantum number. 



1 



ri - V2\ 



s=0 



K(ri,r2) Ps(cos7) , 



(5) 



where 7 is the angle between ri and r2, the Pg are the 
Legendre polynomials and Vs{ri, r2) is given by 



V;(ri,r2) 



for ri < r2 

''2 ^ 

for ri > r2 



(6) 



The expansion in Eq. ([5]) becomes accurate only in the 
limit s 00. However, it is known to be applicable using 
a small Smax as an upper limit of the sum being there- 
fore an expansion parameter. Actually, if the ansatz of 
Eq. ([3|) is used for ip then Smax = 2Zniax holds by what 
^max becomes the decisive expansion parameter. More 
details of the code, which was used in order to calculate 
photon-induced processes [H, \m 
[43] . It is based on an atomic code 
quently applied before, e.g., in calculations of antiproton 
collisions with atomic targets jSlj, ^] . 
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were discussed in 
4711 , which was fre- 



B. Impact parameter approximation 

The collision process is considered in a semi-classical 
way using the impact parameter method (cf., e.g., Ref. 
(48l |) which is believed to be highly accurate for im- 
pact energies E > I keV. The quantum-mechanically 
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treated electron is exposed to the Coulomb potential of 
the molecular nuclei as well as the heavy projectile. The 
latter is assumed to move on a straight classical trajec- 
tory R(i) = h + vt given by the impact parameter b 
and its velocity v while t is the time. The space-fixed 
coordinate system is defined with the x and z axis being 
parallel to b and v, respectively. 

For a fixed Rnuc the collision process can be described 
by the time-dependent Schrodinger equation 

z^vl,(r, R(t)) ^ + V^„t(r, R^)) ^{r, R{t)) , (7) 

where the interaction between the projectile with charge 
Zp and the target electron is expressed by the time- 
dependent interaction potential 

The interaction of the projectile with the nucleus leads 
only to an overall phase which does not change the total 
cross sections. It is therefore not considered in this study. 
The time-dependent scattering wave function 

*(r, R{t)) = J2 Cfc(R(i)) Mr) e'^'-' (9) 

k 

is expanded in the normalized time-independent eigen- 
states as given in Eq. ([3]) and k = A^MII stands for 
the quantum numbers needed to label these states. Sub- 
stitution of Eq. ([3]) into Eq. ([7]) and projection with i/jk 
leads to the usual coupled equations 

for every trajectory R(t), i.e., for every impact parameter 
b and every impact energy E — (1/2) Mpv'^, where Mp is 
the projectile mass. The Cfe(R(i)) depend of course also 
on (the fixed) Rnuc- The differential equations (jTU]) are 
integrated in a finite z-range —50 a.u. < z — vt < 50 a.u. 
with the initial conditions Ck{'R{ti——50/v)) = Ski, i-e., 
the target is initially in the electronic state ipi. 

The probability for a transition into the electronic final 
state ipk at tf = 50/v for a fixed Rnuc = (^nuc, ©, ^) is 
given by 

Pk{b,E;Rnnc,e,<^) = |cfc(&,w,i/;i?„nc,e,$)p. (11) 
In accordance with [41], the transition probability 

Pkib,E)^ J |x.,(i?nuc)i;"(e,<i>)p ^^^^ 

X pk{b, E; Rnuc, 6, $) sinedi?„ucded4> . 

becomes orientation-independent by integration over 
Rnuc- The corresponding cross section 

akiE)^ 27: [ pk{b,E)bdb, (13) 



can then be obtained by integration over b as it is done 
for atomic targets which are spherical symmetric. The 
total cross sections for ionization, 

^ion(£;) = 5] <Tfe(£;), (14) 

f-k>0 

and for excitation of the target, 

«o<«fc<0 

can be obtained by summing up all partial cross sections 
into states k (as given in Eq. (I13p ) with positive energy 
and all Cfe for states with negative energy being larger 
than that of the ground state eo, respectively. 

C. Franck-Condon approximation 

The dependence of the ionization cross section on the 
internuclear distance i?nuc for p + Hg" was examined in 
[4l| for the range 1.5 a.u. < i?nuc < 3 a.u. in which the 
radial distribution |xi^j|^ of the vibrational ground state 
Xoo is non-negligible. It has been shown that the de- 
pendence of the cross sections on i?nuc is approximately 
linear. A similar dependence of the cross sections on 
_Rnuc was also obtained in calculations for p -I- H2 in [38| . 
Under the assumption that is an even function of 

-Rnuc - -Rnuc and cr(i?nuc) is linear in i?nuc around ^nuc 
the Franck-Condon (FC) approximation becomes accu- 
rate as discussed, e.g., in [49|, where -Rnuc = (-Rnuc) is 
the expectation value of i?nuc- Consequently, in [4lj the 
FC results were found to be very close to the exact cross 
sections obtained by an integration over i?nuc like in Eq. 

m- 

In what follows the FC approximation is used, i.e., the 
calculations of the ionization and excitation cross sec- 
tions are performed for i?nuc = 2.05 a.u. which is the 
expectation value for the ground state. 

D. Molecular orientation 

In contrast to atomic targets one set of trajectories 
in which only b is varied is not sufficient for molecu- 
lar targets. In the case of a molecule also different an- 
gular orientations between the trajectory of the projec- 
tile and the molecular axis have to be considered in or- 
der to perform the integration in Eq. p^ . This leads 
to a three-dimensional set of trajectories which can be 
characterized by {6, O, $}. In the space-fixed frame, de- 
fined by b and v, the position vector of the projec- 
tile is R(t) = {Rx, Ry, Rz) — {b,0,vt) and the spher- 
ical coordinates of one molecular nucleus are given by 

(i?nuc/2,e,$). 

The electronic wavefunctions ^ in Eq. ([3]), however, are 
defined in the molecule-fixed frame in which the z' axis is 
chose to be parallel to the internuclear axis (cf. Fig. [1]). 
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(ii) 



Molecule-fixed frame 



(iii)^ 



Figure 1: Sketcli of the molecule-fixed frame. The z' axis is 
parallel to the internuclear axis of the H^. Shown are three 
trajectories of the p and the corresponding impact parameters 
b for the orientations (z), (ii), and (in), which are further 
explained in the text. 



Therfore, also the time-dependent calculations of the col- 
lision process are performed in the latter frame. Therein, 
the position vector of the projectile R = {Rx' , Ry' , Rz') 
can be written as 



Rx 



Rv 



b cos 9 cos ^ — vt sin Q , (16) 
-6sin$, (17) 
R^, = 6 sin6 cos$ + -wi cos6 , (18) 

or be expressed in spherical molecule-fixed coordinates 



R 



Hi, 



arccos 



arctan 



Rz^ 
R 

Rx' 



Rl, = Vb^ + vH^ 



(19) 
(20) 

(21) 



where 4> can take all values within the interval [0, 27r]. 
Note, in Eqs. (|19 [) - (PT|) all three spherical coordinates are 
explicitly dependent on the time as well as on the impact 
parameter b. The angular coordinates 9 and depend 
also on the relative orientation between the trajectory 
and the internuclear axis given by Q and i? = |R| is 
obviously the same in both frames. 



E. Interaction matrix elements 

The matrix elements of the time-dependent interac- 
tion potential induced by the movin g p rojectile are sim- 
ilar to those for atomic targets in [31i|. However, due 



to the loss of the spherical symmetry of the target two 
major differences exist. First, different orientations of 
the molecular axis lead to different interaction matrix el- 
ements. Second, the reduction of symmetries results in 
other (good) quantum numbers {M, IV) and therefore to 
a different scheme of calculating the matrix elements. 

The matrix elements of the interaction potential be- 
tween the two states 4'N'n' ^^'^ V'jvn which are required 
in Eq. pO)l are given by 



M' 



mnt 



M 

Nn 



Vn'w 



5-0 q 



E 



47r 



(2s+l) 



F^Yr(0,0)YKc.) 



(22) 



'max 

~ ^ M'.M / , 



^ M'N'l',MNLs ^M'l'.MLs A^'J; 



i'=i' . ;=;„i„s=sio. 



The three terms V, 
which depend on 
as 



(0) 



M',M' 



V'-'^'^ and 



R and 9, respectively are defined 



(0) 



M'M 



-i(M'-M)4> 



V, 



M'N'l\MNl,s 



M' 
Pn'V 



Vs{r,R)\p'^i) , 



(8) 



M'l'.AILs 



= (-1)*^ ^{21' + l){2l + l) 



(24) 
(25) 
(26) 



{s + {M'-M)y. ' ^ 



I' s I 




/' s I 

-M' (M'-M) M 



where the P™ are the associated Legendre polynomi- 
als. Due to the Wigner-3j symbols in Eq. ([26|) s only 
takes even or odd values depending on whether Sup — 
^max + ^max IS cven Or odd. The lower limit of s is deter- 
mined by Slow = max(|M'-M|, \l'^^^-lnia^\). Note, that 
expression (|23p factorizes into two parts which can be de- 
termined independently, i.e., vlf,\j and a second term 

depending on R and 9 in which different F'-^^ and V^^^ 
are mixed by the summations over Z', / and s. Further- 
more, the behavior of the expressions in Eqs. (|24p - (l^^ 
under exchange of the initial and the final state as well 
as under a simultaneous change of the signs of AI' and 
M 



V, 



i4>) 



M'M 



V, 



M,M' 



V, 



(R) 



M'N'l',MNl,s 
(R) 



v_ 



^-M',-M 

-M'N'l',-MNl,s 
yi^R) 



-M.-M' 



V, 



^ MNl,M' N'V ,s — ^ -MNl,-M'N'l',s ' 
(«) _( -,\( M'-M) Trie) 



M'l',Ml,s 
(6) 



V, 



MLM'V ,s 



{^^l)(M'-M)y_ 



M'l',-MLs 



(27) 
(28) 

(29) 



(6) 

Ml.-M'l',s ' 



can be used to reduce the computational effort. 

While and F^"^) can be determined analytically 
the radial part is integrated numerically using 
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quadrature. Furthermore, the number of different y(^) 
is much larger than those of V'^^^ and y'-'^-' since the for- 
mer depends on all parameters. Therefore, V^-^'^ has to 
be evaluated efficiently. 



F. Symmetries and selection rules 

One aim of the present study is to reduce the com- 
putational effort for p + collisions drastically. This 
permits the description of excitation cross sections but 
it is even more a precondition in view of future calcula- 
tions for the much more demanding system p -I- H2 . As 
can be seen from Eq. (fTT|) a large number of calculations 
are necessary in order to perform an integration over the 
angles Q and $ in Eq. (fT2|) . An alternative approach. 



though approximate, is to use for fixed -Rnuc, b, and E an 
orientationally-averaged transition probability defined by 

1 



Pk = 



Pfe(0,0)+pfc(|,0) , (30) 



in which only the three perpendicular orientations (O, $) 
= (i) (0,0), (m) (7r/2,0), (m) (7r/2,7r/2) are considered 
instead of performing the integration in Eq. (jl2p . 

In this work it was found that the integrated ionization 
cross sections by Sakimoto 4lj — obtained according to 
Eq. — can be reproduced nicely using — according 
to Eq. (pO|) — his results only for the three perpendicular 
orientations (i), (u), and {Hi). The relative difference of 
the ionization cross section obtained by integration and 
by averaging of the probabilities for the three orientations 
is around 1% for £; = 2 keV and 2% for E = 100 keV. 
Therefore, in what follows only the three perpendicular 
orientations (i), (m), and {Hi) are considered, although 
the present method is capable of arbitrary angular ori- 
entations. In Fig. [T] the trajectories for (i), (m), and 
(Hi) are sketched in the molecule-fixed frame. In the fol- 
lowing, the properties of these three trajectories and the 
symmetries of the according interaction matrix elements 
in Eq. (^5)1 are discussed. 

In (i) the molecule- and space-fixed frame coincide re- 
sulting in the same kind of problem as for atomic targets 
(cf., e.g., HH, 113), i.e., cos6' = {vt)/R and cj) = 0. In 
(m) the molecule is oriented parallel to the impact pa- 
rameter which gives equal to or tt for R^' > or 
Rx' < 0, respectively and for 9 one gets cos 9 — b/R. 
Finally, in (Hi) the molecular axis is oriented perpendic- 
ular to the collision plane resulting in a time-dependent 
(j) = arctan(fe/(u t)) while 9 is constant with cos6' — 0. 

For (i) and (m) the azimuthal angle ((> can be consid- 
ered as constant [58*1. As a consequence, the system of 
coupled equations in Eq. (|10p can be transformed in such 
a way that only positive M quantum numbers have to be 
treated explicitly when solving the coupled differential 
equations. Such a transformation was demonstrated in, 
e.g., Ref. [131 ■ Alternatively, the angular part of basis 
states can be described with a combination of spherical 



harmonics ((— l)™!^™ -|- ™) which is solely real as it 
has been done for atomic targets (cf., e.g., Refs. [3ll.[48j). 

In (iii) these simplifications are not possible since 4> is 
time-dependent. Consequently, positive as well as nega- 
tive M quantum numbers have to be considered. How- 
ever, the fact that cos^ = holds can be exploited. As a 
consequence, in the interaction matrix elements of T^'^' 
in Eq. (I^S)) all odd associated Legendre polynomials, i.e., 
s + M' + M being odd, vanish. As a result a selection 
rule only allows for transitions in which the parities of 
the initial and final state differ and the difference of the 
initial and final M is odd or both parities are equal and 
the difference of the M is even. In the case that the 
molecular ion is initially in its ground state only tran- 
sitions among the symmetry subspaces (M, 11) = (0,.g), 
(1,u), (2,ff), (3,m), (4, g),... are allowed. 

Due to the mentioned symmetries for [i) and (m) as 
well as {iii) the eigenstates separate into two sets which 
can be treated independently since they are not coupled 
by the matrix elements of the Coulomb interaction (Eq. 
((23)) ). By this the numerical effort can be reduced by 
nearly a factor 4. Note, although the time propagation 
in {i) and {ii) has only to be performed for either M > 



Table I; Variation of (7ion(0,"I>) and crexc(6, "I"), in units of 
10~^^ cm^, for p-l-Hjwith respect to A'^, the number of B 
splines, while M = 3 and A/ = 2 are kept constant and 
-Rnuc = 2.0 a.u. 



aion(e,$) 



E (keV) 2 



50 



250 



c(e,$) 

50 250 



N 






(e,$) 


= (0,0) 






14 


0.0737 


0.4065 


0.1867 


0.8755 


1.4922 


0.7251 


21 


0.1194 


0.4064 


0.1856 


0.8443 


1.4633 


0.7172 


30 


0.1383 


0.4053 


0.1869 


0.8247 


1.4638 


0.7177 


35 


0.1396 


0.4059 


0.1872 


0.8226 


1.4651 


0.7187 


50 


0.1396 


0.4058 


0.1871 


0.8228 


1.4651 


0.7186 


65 


0.1396 


0.4059 


0.1871 


0.8228 


1.4651 


0.7187 


N 






(e,$) = 


= (^/2, 0) 






14 


0.0340 


0.4574 


0.1980 


0.3165 


2.1410 


1.2231 


21 


0.0447 


0.4453 


0.2023 


0.2994 


2.1068 


1.2137 


30 


0.0476 


0.4433 


0.2031 


0.2973 


2.1074 


1.2127 


35 


0.0480 


0.4439 


0.2032 


0.2980 


2.1072 


1.2127 


50 


0.0480 


0.4438 


0.2031 


0.2979 


2.1073 


1.2127 


65 


0.0480 


0.4438 


0.2031 


0.2979 


2.1073 


1.2127 


N 






(6,$) = 


{■k/2,^/2) 




14 


0.0280 


0.2982 


0.1661 


0.1892 


0.7161 


0.5202 


21 


0.0393 


0.2787 


0.1615 


0.1709 


0.7013 


0.5159 


30 


0.0439 


0.2786 


0.1619 


0.1670 


0.7034 


0.5169 


35 


0.0444 


0.2793 


0.1622 


0.1675 


0.7061 


0.5182 


50 


0.0444 


0.2792 


0.1621 


0.1675 


0.7059 


0.5181 


65 


0.0444 


0.2792 


0.1621 


0.1675 


0.7061 


0.5182 
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or M < the matrix elements in Eq. ((23l) have to be 
computed for negative and positive M which, however, 
differ only in l/t^' and V^'f'\ 



III. RESULTS 

A. Convergence behavior 

In Tables U HlIII the variations of the quantities 
crion(0,<i') and (Toxc(0,'f) with respect to the basis 
set parameters N , M and are presented consider- 
ing the three different orientations (6,$) = (i) (0,0), 
[a) (7r/2,Q), and (m) (7r/2,7r/2)- The crion(e,$) and 
a£xc(6, are defined in accordance with Eq. (14) of Ref. 

m by 

(e,$) = 27r y"pio„(6;e,$)M6, (31) 
(e,$) = 27r y"pcxc(6;e,$)6d&, (32) 

where pion and poxc are the probabilities for ionization 
and excitation, respectively. A fixed internuclear dis- 
tance is used which is chosen to be i?nuc — 2.0 a.u. for the 
convergence study being the equilibrium distance of an 
molecule. Note, that the o'ion(0,^) and croxc(0,'I') 
as given in Eqs. pip and (I32p are no measurable quan- 
tities and are only defined in order to learn more about 
the orientational dependence. 



Table 11: Variation of CTion(6,'l') and (Joxc(6, $), in units of 
10"^*^ cm^, for p-(-Hjwith respect to M while iV = 50 and 
AZ = 3 are kept constant and i?nuc = 2.0 a.u. 







aion(e,$) 


CT" 


oxc(e,'l>) 


E (keV) 2 


50 


250 


2 


50 


250 


M 






(e,$) 


= (0,0) 






1 


0.1371 


0.4118 


0.1547 


0.8046 


1.4858 


0.7184 


2 


0.1353 


0.4045 


0.1839 


0.7936 


1.4404 


0.7062 


3 


0.1350 


0.4004 


0.1889 


0.7916 


1.4367 


0.6998 


4 


0.1354 


0.3999 


0.1894 


0.7923 


1.4366 


0.6982 


M 






(G,$) = 


= (V2,0) 






1 


0.0497 


0.4447 


0.2047 


0.2787 


2.1237 


1.2262 


2 


0.0465 


0.4435 


0.2138 


0.2875 


2.0899 


1.2144 


3 


0.0458 


0.4413 


0.2194 


0.2890 


2.0856 


1.2096 


4 


0.0460 


0.4404 


0.2203 


0.2898 


2.0851 


1.2092 


M 






(9,$) = 


(7r/2,^/2) 




1 


0.0472 


0.2019 


0.1149 


0.1496 


0.7817 


0.4905 


2 


0.0426 


0.2573 


0.1501 


0.1537 


0.6984 


0.4991 


3 


0.0411 


0.2704 


0.1617 


0.1552 


0.6666 


0.4964 


4 


0.0410 


0.2711 


0.1659 


0.1552 


0.6592 


0.4949 



Table U shows the excellent convergence behavior of 
the it(8, $) with respect to the number N of states per 
M, n, and Z, independently of the impact energy, an- 
gular orientation or whether ionization or excitation is 
considered. A relatively small value of iV = 35 yields 
already relative errors |Act(6, $)/cr(8, $)| < 0.1% for aU 
considered crion(6, $) and crcxc(6, <&) of Table HI 

The variation of the (Tion(8,$) and fTcxc(0,$) with 
respect to M in Table [IT] yields relative errors for M — 2> 
which are smaller than 0.5% except for {iii) where these 
are 2.5% for CTion(7r/2, 7r/2) ai E = 250 keV and 1.1% for 
crex(7r/2, 7r/2) at £; = 250 keV. 

A somewhat worse convergence behavior can be ob- 
served in Table IIIII for the variation of A/ = (^max — 
^min)/2 which is the number of different / per symmetry 
subspace (Af, 11) . AZ = 3 gives relative errors less than 
1%, 1.5%, and 2.5% for the orientations (i), [ii), and 
{iii), respectively. Increasing AZ to a value of 4 decreases 
the maximal relative error to 0.6%. In general larger AZ 
are required for lower impact energies. This trend is also 
known from atomic targets [1, [3l| 

In the calculations of the present results the basis set 
parameters are chosen to be {N,M,Al) = (35,3,3) for 
the orientations (i) and (ii) and {N,M,Al) = (35,4,3) 
for (iii). The size of the basis is in general given by 
N X (2M-I-1) X 2AZ. Exploiting the symmetries discussed 
in Sec. Ill Fl a total number of 840 and 945 coupled differ- 



Table III: Variation of (Tion(S,$) and crexc(0, "J"), in units 
of 10~^^ cm'^, for p-l-Hjwith respect to AZ = (Zmax — Zmin)/2 
while A'^ = 50 and M — 3 are kept constant and i?nuc = 2.0 
a.u. 

aion{Q,^) (Tcxc(6,'l>) 

E{keV) 2 50 250 2 50 250 



AZ 






(e,$) 


= (0,0) 






1 


0.2010 


0.3311 


0.1716 


0.7733 


1.5828 


0.7971 


2 


0.1396 


0.4058 


0.1871 


0.8228 


1.4651 


0.7186 


3 


0.1349 


0.3977 


0.1901 


0.7917 


1.4374 


0.6986 


4 


0.1338 


0.3938 


0.1900 


0.7847 


1.4314 


0.6929 


5 


0.1332 


0.3929 


0.1900 


0.7827 


1.4295 


0.6912 


AZ 






(9,$) = 


= (7r/2,0) 






1 


0.0672 


0.3594 


0.1252 


0.3312 


2.0894 


1.0544 


2 


0.0480 


0.4438 


0.2031 


0.2979 


2.1073 


1.2127 


3 


0.0457 


0.4382 


0.2213 


0.2891 


2.0859 


1.2070 


4 


0.0450 


0.4344 


0.2243 


0.2859 


2.0822 


1.2023 


5 


0.0449 


0.4334 


0.2248 


0.2848 


2.0809 


1.2007 


AZ 






(0,$) = 


(^/2,^/2) 




1 


0.0855 


0.4475 


0.1944 


0.2809 


1.1539 


0.7686 


2 


0.0444 


0.2792 


0.1621 


0.1675 


0.7059 


0.5181 


3 


0.0411 


0.2701 


0.1616 


0.1552 


0.6666 


0.4960 


4 


0.0401 


0.2664 


0.1601 


0.1513 


0.6559 


0.4898 


5 


0.0399 


0.2656 


0.1599 


0.1503 


0.6529 


0.4881 
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ential equations have to be solved for (i), (ii) and {in), 
respectively. A further increase of one of the parameters 
N, M, or Al leads for nearly all energies and orienta- 
tions to decreasing (Tion(0, <&) and o-oxc(0, *&)• Therefore, 
the present results obtained with fixed sets of {N, M, Al) 
might be considered as upper bounds to the exact values. 
An energy cutoff of 25 a.u. is used, i.e., only (continuum) 
states V'fc with < 25 a.u. are considered in the expan- 
sion of ^1 in Eq. 



B. Total cross sections 

Calculations for p collisions with are performed 
within a broad energy range of 0.5 keV > i? > 10 MeV. 
The three orientations of the molecular axis in the space- 
fixed frame (e,$) = (0,0), (7r/2,0), and (7r/2,7r/2) are 
considered. Trajectories for these three directions in the 
molecule-fixed frame are sketched in Fig. [T] The FC ap- 
proximation is employed throughout with i?nuc = 2.05 
a.u. An orientationally-averaged transition probability 
Pk (6, E) is gained according to Eq. ((30)) from the results 
for the three orientations. Subsequently, the total ioniza- 
tion and excitation cross sections are obtained as given 
in Eqs. (Unil-dlSD. 

The present data for the total ionization and excitation 
cross sections Cion and ctoxc, respectively, are listed for a 
selection of energies in Table IIVI 



Table IV: Ionization and excitation cross sections for antipro- 
ton collisions with in 10~^® cm^ which are shown in Figs. 
^a.) and[2jb), respectively. The results are given for the mean 
value of the internuclear distance iinuc ~ { Rnuc ) = 2.05 a.u. 
in accordance with the Frank-Condon approximation. The lit- 
erature data for ionization were calculated by Sakimoto [4ll ]. 



E (keV) 


^ ion 




cxc 


1 


0.036 




0.235 


2 


0.077 


0.078 


0.445 


4 


0.147 




0.785 


5 


0.173 




0.909 


8 


0.239 




1.176 


10 


0.268 


0.268 


1.288 


20 


0.348 


0.349 


1.504 


40 


0.383 




1.516 


50 


0.379 


0.380 


1.475 


80 


0.350 




1.373 


100 


0.326 


0.333 


1.271 


200 


0.228 


0.232 


0.940 


400 


0.135 




0.631 


500 


0.112 


0.113 


0.538 


800 


0.073 




0.392 


1000 


0.060 




0.392 


2000 


0.031 




0.189 


4000 


0.016 




0.102 


8000 


0.008 




0.053 



1. Ionization 

The orientationally-averaged FC cross sections for ion- 
ization Uion are shown in Fig. [2jaO and are compared to 
results calculated by Sakimoto [4l| . The present calcula- 
tions reproduce the latter results nearly perfectly. From 
Table IIVI it can be seen that for E < 50 keV the agree- 
ment is better than 1%. For E = 100 keV, 200 keV, and 
500 keV, the difference is of the order of 2.1%, 1.7%, and 
0.9%, respectively. The increased differences between the 
two calculations for the latter three energies might be 
caused by a reduction of the basis size in ^Ij for E > 50 
keV while it is kept the same in the present work. As in- 
tended in [iH this reduction seems to have less influence 
on the results with increasing E but is still visible for 
E < 200 keV. The author of ^4r| also mentions that for 
E ~ 100 keV and (0, <i>)=(0, 0) the relative convergence 
error is largest and not below 2%. 

Due to the good agreement between both calculations 
the following conclusions can be drawn. First, the re- 
sults by Sakimoto [4l[ are confirmed by the use of a sub- 
stantially independent approach. Second, the use of only 
three orthogonal orientations of the molecular axis seems 
to be sufficient for the description of the total ionization 
cross section in p + collisions. Consequently, the 
effort is reduced drastically since a two-dimensional inte- 
gration over the angles Q and $ is not performed. Such 
an integration as given in Eq. p2p requires a sufficient 



number of supporting points, i.e., calculations, in Q and 
$ direction which is of course much larger than three as 
used in the present approach for simple averaging in Eq. 



2. Excitation 

The orientationally-averaged FC cross sections for ex- 
citation ijcxc are shown in Fig. (Hb). To the best of the 
authors' knowledge there are no literature data for ex- 
citation in j? -|- collisions available to compare to. 
Converged excitation cross sections especially for high 
E require an extended range of the impact parameter 
b in comparison to ionization as can be seen in Fig. [31 
In contrast to [4l|, the extention of the b range is well 
feasible with the present approach due to its seemingly 
higher efficiency. In all calculations 30 different b values 
are considered whereas the maximal b increases from 15 
a.u. to 30 a.u. from the lowest to the highest impact en- 
ergies. The spacing between the b values increases with b 
in order to sufficiently resolve the inner region. Fig. [3{b) 
shows for example that a range of 6 < 10 a.u. is not large 
enough for calculating the excitation cross section for an 
impact energy oi E = 125 keV. On the other hand in the 
case of ionization a range of 6 < 5 a.u. is according to 
Fig-ISIa) already sufficient. 
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0.5n-rrTTi I I I 




Energy (keV) 



Figure 2: (Color online) Cross section for p+Hj as a func- 
tion of the antiproton impact energy. Black solid curve with 
plus, present results (orientationally averaged); green circles, 
calculation by Sakimoto 'JT]. (a) Ionization, (b) excitation. 



0.2 I • , • , • , • , r 




Impact parameter b (a.u.) 
0.3| 




Impact parameter b (a.u.) 

Figure 3: (Color online) The transition probability p{b) and 
bp{b) — weighted with b — for p+Hj collisions as a function 
of the impact parameter b for different molecular orientations 
and E = 125 keV. p{b): red dotted curve, (6,$) = (0,0); 
blue solid curve, (7r/2,0); green dashed curve, (7r/2,7r/2). 
bp(b): red circles, (0,0); blue squares, (7r/2,0); green trian- 
gles, (7r/2,7r/2). (a) Ionization, (b) excitation. 



C. Dependence on the molecular orientation 



The shape of CToxc in Fig. ^b) is similar to that of CTion 
in Fig. [2)1 a), although a little less symmetric. The abso- 
lute height of the maximum is, however, about a factor 
4 larger for excitation than for ionization. This factor 
actually is minimal around the maximum and enlarges 
to about 6.5 towards the smallest and largest impact en- 
ergies covered in the present work. The positions of the 
maxima are around 40 keV and 25 keV for ionization and 
excitation, respectively. This is similar to the findings for 
hydrogen atoms but larger than for alkali-metal atom tar- 
gets 3l| . Calculations for p + 1^2 HH in which the 
target was described by an atomic effective one-electron 
model po| yielded maxima for ionization and excitation 
which lie at lower and higher energies, respectively, than 
the present maxima for H^. 



The (Tion(6,$) and o-oxc(6,$) as given in Eqs. ([51]) 
and ([5^ are presented in Figs. IHa) and Hfb), respec- 
tively, as a function of the impact energy together with 
the orientationally-averaged FC cross sections for com- 
parison. The three perpendicular orientations (z), (m), 
and {iii) are considered. 

It can be said that the curves differ considerably for 
different orientations as well as from the orientationally- 
averaged curve for both ionization and excitation. Sim- 
ilarities in the dependence on i?, however, can be found 
for the same orientation between curves for ionization 
and excitation. Thereby, curve (m) and {Hi) share quali- 
tatively a similar behavior while quantitatively the latter 
is for most energies clearly larger than the former. 

For very low energies the curves for {ii) and {iii) seem 
to coincide for ionization {E < 2 keV) as well as for ex- 
citation {E < 1 keV). For these low energies the transi- 
tion probability in Eq. pip is only non- vanishing for very 
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0.5n-rrTTi I I I 




Energy (keV) 



Figure 4: (Color online) Ionization cross section Uion and exci- 
tation cross section CTcxc for p+Hj as a function of the antipro- 
ton impact energy for different molecular orientations. Black 
pluses, orientationally averaged; red circles, (6,$) = (0,0); 
blue squares, (7r/2,0); green triangles, (7r/2,7r/2). (a) Ioniza- 
tion, (b) excitation. 



small b. The differences between the trajectories of the 
p for (ii) and (iii) increase with b but become negligible 
for small b. This can be seen from Eqs. (fTC|) - P^ which 
go for both orientations (m) and {iii) in the limit 6^0 
over to Rx' — —vt, Ry> — 0, and Rz> = 0. Consequently, 
the curves of the transition probabilities in Figs.[2la) and 
[3Jb) for (m) and {iii) merge for b ^ 0. These two curves 
are, on the other hand, most different around b ~ 1 a.u. 
where the trajectories of orientation {iii) encounter the 
position of the nuclei. 

For high energies E > 100 keV, on the other hand, the 
equality 

2aion(0, 0) = aio„(7r/2, 0) + aion(vr/2, Tr/2) = 2f7io„ (33) 

holds with only about 1% deviation. In the case of exci- 
tation another equality, 

croxc(7r/2, 0) = CToxc(0, 0) + croxc(7r/2, 7r/2) = I.Sctcxc , 

(34) 



holds already for E > 50 keV with the same accuracy as 
the one for ionization except for the energies 160 keV < 
E < 250 keV where the deviation is of the order of 5%. 
o'cxc(0,0) shows some structures for energies above the 
maximum in contrast to the orientations (m) and {iii). 

For even higher energies E > 2000 keV all ioniza- 
tion curves seem to coincide. In the case of excitation 
CTcxc(0,0) w CToxc(7r/2,7r/2) holds for E > 4000 keV. The 
behavior at high energies can be explained by the fact 
that the contribution to excitation and ionization from 
distant encounters, i.e., larger b values, increases with 
E. At larger distances from the center the electron is 
exposed to a quasi-central potential and the near-field 
details including the spatial distribution of the nuclei are 
not that much resolved any more. This is especially true 
for the orientations (?) and {iii), as can be nicely seen in 
Fig. [IJb), for which the molecular axis lies in the plane 
perpendicular to b. This is in contrast to {ii) where the 
molecular axis is parallel to b and therefore the minimal 
distance between the antiproton and one of the nuclei 
is smaller. 

The observed dependence on the orientation at high 
energies suggests that ionization and excitation for p + 
can be described in accordance with Eqs. ((33l) and 
(1341) surprisingly well with only one trajectory, i.e., for 
ionization and E > 100 keV by aion — o'ion(0, 0) and for 
excitation and E' > 50 keV by CTcxc — 0.5croxc(7r/2, 0). It 
also means that for £' > 100 keV an appropriate one- 
center model potential might be sufficient in particular 
for the description of the ionization process. For these 
high energies satisfactory results were obtained in calcu- 
lations for p + II2 collisions [H, |5l[ using a one-center 
one-electron model for the description of the H2 molecule 
4C|. The calculations reproduced the experimental data 
lll33j for £^ > 90 keV. For lower energies the mentioned p 
-I- H2 results resemble qualitatively those for the orienta- 
tion {i) which separate from the orientationally-averaged 
curves in Fig. 2] for E < 100. Note, only one kind of tra- 
jectory is possible with the employed H2 model potential 
due to its atomic (spherical-symmetric) character. This 
kind of trajectory is practically the same as the one for 
the orientation (i) in the molecule-fixed frame. 

Finally, it is interesting to note that the knowledge 
of the results for maximally three perpendicular orienta- 
tions appears to be sufficient to accurately reproduce the 
total angular-integrated ionization cross section although 
the three curves differ considerably. 



D. Production of H+ in p + and p + H2 

1. P + Ht 

In collisions of p with three main mechanisms lead 
to the production of H"*". First, ionization of the target 

P + H+ ^p + H++H+ + e-, (35) 
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second, electronic excitation followed eventually by dis- 
sociation 

p + H+^p + H + H+, (36) 

third, direct dissociation which is not included in the 
present approach. At very low energies the replace- 
ment of the electron by the p may also become impor- 
tant. However, these energies lie beyond the scope of the 
present work [s^, 113 ■ A_s can be seen from Eqs. ([55)1 
and (|36p two are produced in the ionization process 
while it is only one for excitation. The cross section for 
H+ production due to ionization and excitation is there- 
fore given by the sum 2crion + fcxc- It follows from the 
present results (cf. Table [TVl and Fig. [2]) that the absolute 
contribution to this sum in p + collisions originates, 
however, only from a third to a quarter from ionization 
depending on whether the impact energy is close to or 
away from the maximum, respectively. 

2. p + H2 

While no experimental data are available for p + 
collisions measurements of the ionization and the H+ pro- 
duction cross sections have been performed for p + H2 
[l|. The present results obtained for p + collisions 
shall be used to learn more about the different produc- 
tion mechanisms of the measured H"'' cross section for p 
+ H2 which has not been explained theoretically so far. 
The analysis is done by applying the lEV as introduced 
in Sec. HI It was used by Wehrman et al. j22j for the de- 
scription of double ionization vap + He collisions. In the 
lEV double ionization is considered as a two-step process 
and the cross section is obtained by using the product of 
transition probabilities from (effective) one-electron cal- 
culations only. 

For the description of the first step the single-electron 
ionization probabilities for p -I- H2 are calculated 
as explained in [38j . Therein, a simple one-center one- 
electron model potential for the description of the H2 
target [40j was used which reproduces experimental ion- 
ization and excitation data well for £' > 90 keV. 

In the second step p + HJ collisions can contribute in 
two ways to the cross section for H+ production (in con- 
trast to p + He^ in double ionization of helium). They 
are given in Eqs. ([55]) and ([M)) each having the proba- 

. . H+ H+ 
bility Pj^Jj and pcxc, respectively. Therefore, both cross 

sections, (Tdi for double ionization and for ionization 

followed by excitation of an H2 target, 

^di = 2^1 pll{b)pfi{h)hdb, (37) 

(Ji, = 2 7T J pllib) pfi{b) b db , (38) 

are considered in accordance with the lEV. 

Furthermore, all doubly-exited electronic states of H2 
are in principle autoionizing. But it is also possible that 




Energy (keV) 

Figure 5: (Color online) Cross sections leading to the produc- 
tion of H"*" in collisions with H2 targets as a function of the 
projectile energy, p impact. Present results, p: red stars, dou- 
ble ionization adi; blue triangles, ionization and excitation uie; 
green plus, summed production crjj+ =2(Tdi-l-o"ie. p: violet 
X, H"*" production a-^+ . Experimental data: black squares, 
total H"*" production for P-I-H2 black circles, dissociative 
ionization for P+H2 [s^. The inset shows a high-energy cut- 
out of the same data. 



the doubly-excited H2 dissociates before an electron is 
ejected. The description of this channel is clearly very 
subtle and has been studied in detail with a considerable 
effort for the excitation of H2 by photons [13, UHl ■ The 
double-excitation cross section (Tdc for p + TI2 collisions, 

dde = 2^ Jpf;^ib)pf^^{b)bdb, (39) 

is approximated using the independent particle model 
[5^ . The single-excitation probabilities p^^^ for H2 orig- 
inate form the same calculations as the p^^^ used in Eqs. 
dSZl) and 

In Fig. [5] the results for and crjc are presented to- 
gether with the sum of these cross sections 

crH+ = 2o-di + Cio , (40) 

in which CTdi is counted twice since two H+ are produced 
in the double ionization of H2 . Also given are the mea- 
sured cross sections for the total H"*" production in p -|- 
H2 by Hvelplund et al. [H and for dissociative ionization 
for p -f H2 by Shah and Gilbody [36*], where the latter 
includes (Xdi, (Tic and transfer ionization. 

For all energies considered in Fig.[5]crdi lies below (Tie by 
about a factor ^ 2. Therefore, both mechanisms in the 
H+ production as given in Eqs. (j35p and (|36p contribute 
with a comparable amount of H"*" in (7fj+ . Although the 
measured data for p -|- H2 have a similar slope than the 
present (Th+ the latter is shifted down relative to the 
experimental results by between w 10"^^ cm^ for the 
lowest and « 6 x 10~^* cm^ for the highest energies in 
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Fig. [5l The cross section for double excitation (Tdc which 
is not included in (Tjj+ is of the order of approximately 
10% of (Th+ in the whole energy range. This means that 
the three mechanisms for the production of H"*" described 
with the employed models are not sufhcient in order to 
reproduce the experimental data of Hvelplund et al. [H . 

Note, the curves given in Fig. [5] are calculated with 
the FC data presented before which implies that the in- 
ternuclear distances of the is set to i?nuc = 2.05 a.u. 
Under the assumption that the time between the first 
and the second step in the lEV is too short to allow for 
a change of the internuclear distance the p + colli- 
sions are also calculated for i?nuc = 1-4478 a.u. which is 
the expectation value of H2. This leads to an increase 
of the binding energy and therefore to a decrease of the 
ionization and excitation cross sections p + . Con- 
sequently, the results for i?nuc = 1-4478 a.u., which are 
not shown in Fig- [51 become even smaller and reproduce 
those for i?nuc = 2-05 a.u. from about 80% for E — 50keV 
to 90% for E = 1000 keV. 

The experimental data for dissociative ionization in p 
+ H2 collisions by Shah and Gilbody also have a similar 
slope as the present results. Note, in contrast to the 
measured p data, their absolute values are comparable 
with those of the present (7^+ for P- 

Additionally, (7^+ results for p collisions with H2 are 
calculated and shown in Fig- [S] They are obtained ex- 
actly in the same way as described for p impact only that 
the projectile charge Zp in the interaction potential (cf. 
Eq- (O) is set to -1-1 instead of -1 for p. Although the 
present approach does not distinguish between ionization 
and electron capture by the proton, the p results are still 
meaningful for high energies since the cross section for 
electron capture for p -I- H2 is negligible for E > 200 keV 
[37I [57} . In general, the present data for p and p impact 
are very similar and practically the same for E > 400 kcV 
both being close to the experimental proton results. This 
means that an obvious difference of the production 
between p and p impacts for high energies as measured 
experimentally and suggested by the double-ionization 
cross sections for He targets cannot be reproduced by 
the present study. 

Within the employed two-step model it might be cru- 
cial to consider an orientational dependence also for 
the p^^^ as is done in the present method. That way, 
the probability pj^^ can be multiplied first for each ori- 
entation individually with the probabilities for Hj, as 
in the Eqs. (|57)) and (|38p and being only afterwards 
orientationally-averaged. In the case that the depen- 
dence on E of the p^^^ for the three orientations is similar 
to that of the target as shown in Fig. d] this might 
load to a sizable effect on tTdi and trie and therefore also 
on crH+ . 

However, an advanced treatment of the p -I- H2 collision 
would be eligible which includes at least a two-electron 
description of the target in contrast to what has been 



done in [38| in order to resolve the discrepancy which 
appears between experiment and theory. 

IV. SUMMARY AND CONCLUSION 

Time-dependent close-coupling calculations using a 
spectral expansion are performed in order to determine 
ionization and excitation cross section for p + HJ colli- 
sions in a broad energy range from 0.5 keV to 10 MeV. 

For the description of the target the Born- 
Oppcnhcimcr approximation and a one-center expansion 
of the molecular potential are used. The collision pro- 
cess is treated within the impact parameter method as 
well as the Franck-Condon approximation. The transi- 
tion probabilities are obtained by averaging the results 
for only three perpendicular orientations of the molecular 
axis with respect to the p trajectory instead of performing 
an integration over all orientations by what the numer- 
ical effort is reduced drastically. The use of symmetries 
for these three orientations leads to selection rules which 
further reduce the effort by about a factor 4. 

Extensive convergence studies assure that the final re- 
sults do not depend on the expansion parameters. The 
present ionization cross section reproduces nicely the re- 
sults calculated by Sakimoto '4l'| showing that the use of 
three perpendicular molecular orientations is sufficient. 
An extension of the range of the impact parameter b 
in comparison with [4l| is possible due to the less de- 
manding calculations. The larger b range allows for the 
determination of the excitation cross section. In gen- 
eral, the contribution to the total cross sections for a 
fixed impact energy differs considerably for the three per- 
pendicular orientations. However, for high energies the 
cross sections for ionization {E > 100 keV) and excita- 
tion {E > 50 keV) can already be generated with only 
one molecular orientation. 

The cross sections for double ionization and ionization 
followed by excitation are studied as well as their con- 
tribution to the H+ production in p -f H2 collisions. A 
sequential two-step model also referred to as independent 
event model is employed to extract these cross sections 
only from (effective) single-electron transition probabili- 
ties. The transition probabilities for Hj targets are taken 
from the present study while those for single ionization of 
H2 are calculated according to [3l|. The present results 
do not match and are smaller than the experimental data 
[l| from 50 keV to 2 MeV. The present data are, how- 
ever, very similar in the case that p -|- H2 collisions (36j 
are considered. This motivates further experimental and 
theoretical work on p -I- H2 collisions. 
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